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Abstract. The azimuthal distributions of photons and charged particles with respect to the event plane are investigated 
as a function of centrality in ^"**Pb -I- ^"**Pb collisions at 158 -A GeV/c in the WA98 experiment at the CERN SPS. The 
anisotropy of the azimuthal distributions is characterized using a Fourier analysis. For both the photon and charged 
particle distributions the first two Fourier coefficients are observed to decrease with increasing centrality. The observed 
anisotropics of the photon distributions compare well with the expectations from the charged particle measurements for 
all centralities. 



PACS. 25.75.Dw Particle production, azimuthal anisotropy, flow 
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1 Introduction 

Non-isotropic emission of particles with respect to the reaction 
plane, as first observed at the Bevalac 1 1 1, provides evidence 
for collective flow in high energy heavy ion collisions. Flow, 
or anisotropic particle emission, has been observed for a large 
variety of interacting systems from incident energies of a few A 
GeV/c at the Bevalac (SIS) and AGS to much greater energies 
at the SPS and RHIC II2I3I4I5I6I7I8I9I10I11I12I13I . 

Anisotropic flow manifests itself as asymmetries in the az- 
imuthal distribution of particles and can be reproduced in theo- 
retical models with different underlying assumptions. One sce- 
nario is incorporated in transport models where the particles 
have a mean free path comparable to the system size 1141151161 . 
The models can describe the observed flow up to AGS energies. 
The other scenario applies when the mean free path is much 
smaller than the system size which allows the description of the 
equilibrated system in terms of macroscopic quantities 1 17 18 1. 
Hydrodynamic models are able to describe the qualitative fea- 
tures of the observed flow 1 19 1 for py below about 3 GeV/c. 

The initial asymmetry in the overlap zone of the colliding 
nuclei translates into unequal pressure gradients in different di- 
rections that leads to an elliptic final state momentum distri- 
bution of the particles 1 17 1, causing an elliptic pattern of flow. 
The elliptic flow is therefore expected to be sensitive to the sys- 
tem evolution at the time of maximum compression |20 1 and is 
shown to be sensitive to the equation of state of the compressed 
nuclear matter The variation of asymmetry with centrality en- 
ables to relate the observed flow to the geometry of the overlap 
region I21I16I . One would then expect a scaling of the data 
from AGS to SPS and RHIC provided the physics of elliptic 
flow remains the same |22|. In the case that there is a phase 
transition from hadronic matter to a quark gluon plasma, it is 
expected that the reflection of this transition in the equation of 
state of the dense nuclear matter would result in changes in the 
pressure gradients which would then be reflected in changes in 
the particle flow pattern. 

The first evidence of azimuthal anisotropy at SPS energies 
was observed in the distribution of photons from Sh-Au colli- 
sions at 200 -A GeV measured in the preshower photon mul- 
tiplicity detector of the WA93 experiment at CERN 1 8 1. Since 
almost 90% of photons produced in ultra-relativistic nuclear 
collisions originate from the decay of tt^'s, the anisotropy of 
the observed photon distributions should reflect the anisotropy 
of the 7r° production followed by the effects of decay of the 
TT^'s. Methods have been proposed to estimate the anisotropy 
of the neutral pion emission by measuring the anisotropy of 
photons |23|. The decay introduces non-flow correlations be- 
tween the photon pairs due to four-momentum conservation 
and may dilute the correlations between the tt'^'s and the event 
plane. Determination of the effect of decay enables the deduc- 
tion of the anisotropy of the neutral pions. The photon aniso- 
tropy measurement thus complements the study of the aniso- 
tropy of charged particle distributions. 

In the present work we report results from the WA98 ex- 
periment on the centrality dependence of the anisotropy coef- 
ficients extracted from measurement of the azimuthal distribu- 
tions with respect to the event plane of photons and charged 
particles in the same pseudorapidity interval. Preliminary re- 
sults on the anisotropy of photon emission in Pb-nPb collisions 



have been reported earlier 124125 1 . The paper is organized in 
the following manner: Section 2 describes the experimental 
setup and data selection. The analysis technique and the re- 
sults on the centrality dependence of the azimuthal anisotropy 
of photons and charged particles are discussed in Section 3. 
Section 4 summarises our investigations. 



2 WA98 Experiment and Data Selection 

The WA98 experiment at CERN 1261 placed emphasis on si- 
multaneous detection of hadrons and photons. The experimen- 
tal setup consisted of large acceptance hadron and photon spec- 
trometers, detectors for photon and charged particle multiplic- 
ity measurements, and calorimeters for transverse and forward 
energy measurements. The experiment recorded data with 158-A 
GeV Pb beams from the CERN SPS in 1994, 1995, and 1996. 
The results presented here are from a portion of the Pb run 
in 1996 during which the magnet (GOLIATH) was turned off. 
The analysis presented here used data recorded with the pho- 
ton multiplicity detector (PMD) and the silicon pad multiplicity 
detector (SPMD). The data from the mid-rapidity calorimeter 
(MIRAC) was used to characterize events on the basis of cen- 
trality of the collision. 

The circular Silicon Pad Multiplicity Detector (SPMD), used 
for measurement of the charged particle multiplicity, was lo- 
cated 32.8 cm from the target. It had full azimuthal coverage in 
the region 2.35 < ?/ < 3.75 (beam rapidity ?;hea,„ — 5.81). The 
detector had four overlapping quadrants, each fabricated from 
a single 300 /iw thick silicon wafer. The active area of each 
quadrant was divided into 1012 pads forming 46 azimuthal 
wedges and 22 radial bins with pad size increasing with ra- 
dius to provide a uniform pseudo-rapidity coverage. The in- 
trinsic efficiency of the detector was better than 99%. During 
the datataking, 95% of the pads worked properly. The SPMD 
was nearly transparent to high energy photons since only about 
0.2% are expected to interact in the silicon. Multiple hits of 
charged particles on a single pad were treated as a single hit for 
the present analysis, and are counted as Nhus- The detector is 
sensitive to all charged particles. The energy cutoff appears as a 
low noise threshold. Details of the characteristics of the SPMD 
can be found in Refs. 1 27 .28J . 

The photon multiplicity was measured using the preshower 
photon multiplicity detector (PMD) located at a distance of 
21.5 meters from the target. The detector consisted of 3 radia- 
tion length (Xq) thick lead converter plates placed in front of 
an array of square scintillator pads of four different sizes that 
varied from 15 mmx 15 mm to 25 mmx25 mm, placed in 28 
box modules. Each box module had a matrix of 38x50 pads 
which were read out using one image intensifier + CCD cam- 
era. Details of the design and characteristics of the PMD may 
be found in Ref. |2930|. 

The clusters of hit pads with a signal above a hadron rejec- 
tion threshold were identified as photon-like. The present anal- 
ysis has been performed with the photon-like clusters, which 
are referred to as photons for brevity. Detailed simulations showed 
that the photon counting efficiencies for the central to periph- 
eral cases varied from 68% to 73%. The purity of the photon 
sample in the two cases varied from 65% to 54%. Most of the 
contaminants of the photon sample are charged particles which 
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Table 1. Centrality selections used in the present analysis based on 
the measured total transverse energy. The corresponding fraction of 
the minimum bias cross section, number of participants, and the aver- 
age photon and charged particle multiplicities measured in the pseudo- 
rapidity interval 3.25 < r) < 3.75 are given for each centrality selec- 
tion. 
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< Nhits > 


40.0- 89.9 
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41.1 


34.3 
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150.5 


298.6-347.6 


5-10 


300.4 


190.1 


177.1 


>347.6 


0-5 


353.4 


222.9 


210.3 



deposit enough energy to fall above the hadron rejection thresh- 
old. The hadron rejection threshold is taken as three times the 
energy deposited by a minimum ionizing particle. For photons 
this leads to a low p^ threshold of 30 MeV/c. 

The transverse energy was measured with the MIRAC calo- 
rimeter [31] located at 24.7 meters downstream from the tar- 
get. The MIRAC was used to measure the total transverse en- 
ergy by measurement of both the transverse electromagnetic 
(E^) and hadronic (E!^"-'^) energies in the pseudorapidity in- 
terval 3.5 < r] < 5.5. The measured total transverse energy, 
Et, provides a measure of the centrality of the reaction. Events 
with large E^ correspond to very central reactions with small 
impact parameter, and vice versa. 

The minimum bias E^ distribution has been divided into 
different fractions of the minimum bias cross section corre- 
sponding to different centrality bins |30 1. The most central se- 
lection corresponds to that 5% of the minimum bias cross sec- 
tion auB with largest measured Et. A total of about 0.25 Mil- 
lion events have been analysed. The minimum number of events 
in any centrality selection is 15K and the maximum is 70K. Ta- 
ble^shows the percentage cross section and the corresponding 
number of participants for each centrality bin. The results pre- 
sented here use only the data for the pseudorapidity region of 
common coverage of the PMD and SPMD (3.25 < tj < 3.75) 
where both detectors have full azimuthal coverage. The aver- 
age measured photon and charged particle multiplicities for this 
region of acceptance are also quoted for each centrality in Ta- 
bled 



3 Analysis 

The anisotropy of the azimuthal distribution of particle emis- 
sion with respect to the reaction plane (or event plane) is char- 
acterized by the coefficients of the Fourier expansion of the az- 
imuthal distribution |32|. The first and the second coefficients 
are measures of the directed and elliptic flow when the expan- 
sion is made about the reaction plane, the plane defined by the 
beam direction and the impact parameter. This may be written 
as 



where cj) is the azimuthal angle of the measured particle and 
tpR denotes the azimuthal orientation of the reaction plane. The 
reaction plane can be most accurately determined in an exper- 
iment that measures the (transverse) momenta of the target or 
projectile fragments. Both t^i and V2 can take positive or neg- 
ative values. By convention, positive (negative) values of vi in 
Eq. [2 denote flow (anti-flow) in the direction of the deflected 
projectile fragment, and positive (negative) values of V2 indi- 
cate in-plane (out-of-plane) flow. 

Though the most accurate determination of the reaction plane 
requires the measurement of target (or projectile) fragments, 
most experiments assume that the measurement of any parti- 
cle type, in any kinematic window enables a determination of 
the reaction plane. We wish to distinguish between the plane 
determined by projectile or target fragments and the plane de- 
termined by any other particle type, and throughout this article 
refer to the latter as the event plane. Obtaining the values of co- 
efficients after projecting azimuthal angles on the event plane 
determined from the same set (after removing auto correla- 
tions) maximises the values of anisotropy coefficients, and may 
include non-flow correlations. The values obtained are neces- 
sarily positive, and have been shown as such in the present 
work. The coefficients determined by projecting on an event 
plane from any other set of particles are expected to be smaller, 
and can also have negative values. The difference in the values 
determined using different sets of event planes have been in- 
cluded in systematic error by various experiments L10.ll.13l . 



3.1 Method 

In the present analysis, the azimuthal distributions of particles 
for any particle species in any pseudorapidity window is ex- 
panded as a Fourier series where the coefficients of expansion 
determine the shape of the event. Retaining terms up to sec- 
ond order coefficient in the expansion, the shape can be char- 
acterized by an ellipse for small values of the coefficients. The 
direction of the centroid and the major axis of the ellipse are 
determined from the azimuthal distributions of the particles. 
These directions, along with the beam direction, define the first 
order and the second order event plane respectively, and are 
obtained as 1 32 1 



^L 



1 

m 



tan 



Swj cos m(j)i 



(2) 



2Tr dN 



= 1 + 2wi cos(0 — ipn) + 2w2cos2(0 — ipnXl) 



where m = 1 or 2 for the first and the second order, respectively. 
The (pi are the azimuthal angles of the emitted particles with re- 
spect to a fixed laboratory direction and the Wi are the weight 
factors. For the azimuthal distribution of the particle yield, as 
in Eq.^ the weight factors are set equal to one. In reality, due 
to finite particle multiplicities tp'i and ^■'2 fluctuate about the 
actual event planes that represent the direction of the centroid 
and the direction of the major axis of the elliptic shape. To the 
extent that initial state nuclear densities are spherically sym- 
metric and the density fluctuations are negligible, the initial 
nucleon density in the overlap region is symmetric about the 
impact parameter or reaction plane and so it is expected that 
the two event plane angles are either the same or perpendicular 
to the reaction plane. 



4 M.M. Aggarwal et al.: Azimuthal Anisotropy of Photon and Charged Particle Emission in ^"^^Pb + ^"^Pb Collisions at 158 -A GeV/c 



The anisotropy, or Fourier coefficients of order n, can be 
determined from the azimuthal distribution of the particles with 
respect to the event plane angle of order m, provided n is an 
integral multiple of m, by fitting to the following equation l32l 



dN 



d{ct^-^'J 



«i + E,° 



2vL 



^'J (3) 



u,'j„j is a measure of the offset of the centroid of the distribution 
when n ■ m= \ and is a measure of the difference between the 
major and the minor axes of the ellipse when n ■ m = 2. The 
actual coefficients are obtained from the observed coefficients 
u^j„ as described later. 

Since the event planes do not depend on the geometrical 
setup of the experiment, the distribution of the event plane an- 
gles determined for a large number of events is expected to be 
uniformly distributed in laboratory angle. Any non-uniformity 
in the acceptance of the detectors over the full azimuth will be 
reflected in a non-uniform distribution of event plane angles. 
Any non-uniformity that is identifiable within a large subset of 
events can be coiTected for by appropriate coiTection methods. 
The method employed in the present work is summarized in the 
following. 



3.2 Detector Acceptance Correction 

The coiTected event plane angle is obtained by shifting the ob- 
served event plane angles tp'^ by Aip'^^ 1321 where the latter is 
written as 



^ 2 
^V'm = y2 (- (sin(nTO V™)) cos(nm '0^„) + 

n—l 

(cos(nm?A;„))sin(nmf/'™)) (4) 

where N = 4/m is sufficient to flatten the raw V' distribution. 
The angular brackets denote an average over all events and are 
obtained from the raw distribution of the m*^ order event plane 
distribution. 

The distribution of the first and the second order event plane 
angles, coiTected for acceptance, is shown in Fig. [2 for the 
charged particle hits in the SPMD and the photon hits in the 
PMD. 



3.3 Event Plane Resolution Correction 

The average deviation of the estimated event plane from the 
true event plane due to multiplicity fluctuations can be deter- 
mined experimentally and is termed as the resolution connec- 
tion factor (RCF). Experimentally, RCF is obtained using the 
subevent method described in II32I . The particles in each event 
are sorted in ascending order of pseudorapidity. Each event is 
divided into two equal multiplicity subevents separated by a 
pseudorapidity interval. The two subevents are separated by 
one pad for the SPMD and by At] = 0.05 for the PMD. The ac- 
tual location of the pad in SPMD and the Arj interval in PMD 
is aUowed to vary in the region 3.25 to 3.75 to ensure equal 



multiplicity for the two subevents. The event plane angle ip'^ is 
determined for each subevent. ' 

This enables determination of a parameter Xm directly from 
the experimental data using the fraction of events where the 
correlation of the planes of the subevents is greater than tt/2 
Il32l33l: 



NeyenUm\^P'^ - <| > 7r/2) _ e"^ 



N, 



(5) 



total 



where Ntotai denotes the total number of events, -ip'^, ip'^ are 
the observed event plane angles of the two subevents (labeled 
a and b) and the numerator on the left denotes the number of 
events having the angle between subevents greater than Tr/2m. 
The parameter Xm so obtained is then used to determine RCF„m 
= {cos{nm{4>'m ~ V'm"^)))' where V'm"'^ is the true direction of 
the event plane, and the average is over all events. The RCF 
can be determined from Xm by the following relation from ref- 
erence 1321 . 



/ii^(X™/4) + /^(X™/4)' 



(6) 



where li, are the modified Bessel functions of order v. 

The errors on the RCF values have been obtained by con- 
sidering that the eiTor on N^vents is statistical. The new values 
of Xm are then calculated for values of Nevents ± V^events 
and used in Eq. |6lto calculate new values of RCF. The change 
in the RCF values gives the statistical error on the RCF determi- 
nation. In general, the eiTors determined in this way are asym- 
metric. Symmetric errors have been quoted using the larger of 
the two asymmetric errors. 



3.4 Anisotropy Coefficients 

The anisotropy coefficients are obtained by filling up parti- 
cle azimuthal distributions of one subevent with respect to the 
event plane of the other subevent where the subevent division 
is described in the previous section. The anisotropy coefficients 
have been determined by three methods which differ in detail 
and provide a consistency check. In this analysis the Fourier 
coefficients prior to event plane resolution correction u'j„j are 
extracted for the case with event plane order equal to the order 
of the extracted Fourier coefficient, i.e. u'j„j = ujj„ which we 
will denote by u'j . 

In the first method we determine v'^ = (cosn((/)° — ipn)) 
and v'J^ = {cosn{(f>^ — ip'n)) where 0" represent the azimuthal 
angles of particles in subevent a and ijj'^ is the event plane an- 
gle determined using particles in subevent b. The averages are 
computed over all particles over all events. In the absence of 
non-flow coiTelations, u„ can be determined using 



{cosn{'p'^ - tp'^)) 



(7) 



' The distribution of the corrected event plane angles for these 
subevents is observed to be flat, for both orders, for photons and for 
charged particles. 
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This determines the magnitude of the coefficients and is neces- 
sarily positive. 

The distributions with respect to the event plane are shown 
in Fig. 121 for both photon and charged particles, for both or- 
ders for the centrality selection corresponding to 10-20% of 
the cross section. In the second method the distributions have 
been fitted to Eq.lsjto determine v[i. The corresponding distri- 
butions with respect to the second order event plane are shown 
in Fig. 12] and have been fitted to Eq. |3]to determine W22- For 
both orders, the fits have been made keeping terms up to values 
of n • m = 2 in the summation in Eq.|3l 

The w^„„ are the values determined with respect to the es- 
timated event plane and must be corrected for the event plane 
resolution 1321 to obtain the actual anisotropy values 



RCFn 



(8) 



The factor \/2 arises because the particle distributions have 
been obtained with respect to the event plane of a subevent 
with half of the total event multiplicity. The subevent resolu- 
tion, ^J {cos{nm{%ljf-^ — "0^))), is averaged over all events and 
is in good agreement with RCF„m/V2. 

In the third method, the values Vnn have been obtained di- 
rectly by the subevent method from Xm of Eq.|5]and the fluc- 
tuation in the average multiplicity M of the full events in that 
centrality bin. 



Xn 



(9) 



The different methods yield consistent results for both the 
first order and the second order anisotropy coefficients and the 
difference in the values is included in the systematic error The 
equivalence of the first two methods arises due to the equiv- 
alence of the geometric mean (method 1) and the arithmetic 
mean (method 2) of the resolution uncorrected values v'^ of the 
two subevents and due to the l/VAf dependence of the fluctu- 
ation in event plane determination. The latter contributes to the 
equivalence of the third method with the first two. The three 
methods may yield different values in the presence of non-flow 
correlations where the differences will be governed by the na- 
ture and strength of these correlations. The most probable val- 
ues are determined using method 1 above. The following sec- 
tions discuss the systematic effects that distort the measured 
anisotropics and the centrality dependence of the anisotropics 
for charged particles and for photons. 



4 Anisotropy in Charged Particles 

The measured anisotropy is expected to be less than the ac- 
tual anisotropy due to the finite granularity of the detector. The 
measured values of anisotropy may also be affected by the im- 
precision in the vertex position due to the finite spread of the 
beam. These effects are particularly relevant for the charged 
particle distribution measurement due to the relatively coarse 
segmentation and close proximity to the target of the SPMD 
detector. The effects are estimated using simulations. 



4.1 Granularity 

The finite granularity of the SPMD detector causes a dilution 
of the anisotropy of the azimuthal distribution primarily due 
to efficiency losses from multiple hits. The quantitative effect 
of the finite granularity on the measured anisotropy has been 
estimated and is briefly discussed in the following. This has 
been corroborated by simulations and is described in greater 
detail in 1 34 1 . 

Defining the mean occupancy as the ratio of the number of 
particles incident on the detector to the number of active cells 



N, 



M 



part 



Nc. 



one can show that 



/i = ln(l + 



N,„ 



Nu 



(10) 



(11) 



cupied cells 1351 . One can further show that 



where Nocc and Nunocc are the number of occupied and unoc- 



(12) 



N, 



part 






where Nhn is the number of ocupied cells. 

Since the intrinsic occupancy of cells increases with the in- 
crease in the number of incident particles, the occupancy will 
have the same azimuthal dependence as the incident particles. 
Substituting this for the occupancy in Eq.^] gives us the az- 
imuthal dependence of the hits. Expressing the resulting equa- 
tion as a Fourier series and collecting coefficients of cos ncf) en- 
ables a determination of the ratio t;^***/t;„, which to first order 
can be approximated by 



,,hits 



1 



ln{l — x) 



(13) 



where x ■ 



Nh 



(= 



A^o 



-). 



The anisotropy in the distribution of charged particles can 
be obtained using the measured value of anisotropy in the dis- 
tribution of hits. 

The azimuthal distributions of the charged particles were 
generated with different initial anisotropics with multiplicities 
corresponding to the measured results. The charged particle 
hits were sorted into the SPMD bins (2° (f>-hm) assuming a 
94% detection efficiency but taking into account hit losses due 
to multiple hits in a single SPMD pad. The resulting azimuthal 
distributions were then analyzed to determine the anisotropy 
coefficients using the method detailed in Sec 3.4 above. Fig.|3l 
shows the results of the simulation for both orders, where the 
ratio of the estimated anisotropy from the hit distribution to the 
initial anisotropy is shown for varying hit multiplicities, far be- 
yond the range of measured values in the SPMD. This is done 
for different values of initial anisotropy. The correction factor 
as given by Eq.^jis shown as a solid line in the figure. One 
observes that the simulation results corroborate the results ob- 
tained above. The simulation results show that the extracted 
anisotropy is systematically lower than the initial anisotropy. 
Part of this loss occurs directly due to efficiency losses from 
multiple hits, and contributes both to the anisotropy and the 
event plane resolution. 
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4.2 Shift in Vertex and Beam Spread 

The finite beam size caused an imprecision in the assumed ver- 
tex by up to one mm in the WA98 experiment. A small shift 
in the vertex position does not affect the azimuthal distribution 
in the fine granularity PMD, situated at 21.5 meters from the 
vertex. However, it can produce an apparent anisotropy in the 
hit distribution in the SPMD situated at 0.328 meters from the 
target. 

The effect of vertex shifts due to the beam spread was also 
investigated by simulation. The vertex position was generated 
according to a two dimensional Gaussian with width a. Parti- 
cles were simulated to originate from this vertex position with 
a realistic rj and px distribution and an azimuthally symmet- 
ric (f) distribution. The particles were projected onto the SPMD 
plane, and their hit positions recorded according to the granu- 
larity and nominal location of the SPMD detector The recorded 
positions corresponded to values of 77 and which differed 
from the generated values due to the shifted position of the ver- 
tex. The simulated distributions were then analyzed in the same 
manner as the experimental data. This was repeated for differ- 
ent values of a and the values of Vn were obtained for each 
sample generated. The maximum possible shift has been de- 
duced by assuming an azimuthally symmetric distribution for 
the most central class, and assigning the granularity corrected 
observed value for first order anisotropy to the shift in vertex. 
This corresponded to a maximum width of the Gaussian distri- 
bution due to beam spread of 0.3 mm. 

A shift in the vertex position due to beam spread, or time 
variation during the 2.5 s SPS spill, produces an azimuthal dis- 
tribution which has a shifted centroid. If the shift occurs on 
an event-by-event basis it cannot be corrected for by the ac- 
ceptance correction methods discussed above, since they can 
only correct for average effects. To investigate the possibility 
of systematic shifts correlated with time during the SPS spill, 
the SPMD charged particle azimuthal distributions were ana- 
lyzed for different times during the spill. No significant varia- 
tions with time during the spill were observed. 

The second order anisotropy is obtained from the fit to the 
elliptical shape of the measured hit distribution. A vertex shift 
does not affect the elliptical shape, as verified by simulations. 



4.3 Results 

The finite granularity requires a correction which is obtained 
using the measured values of u„ from the hit distributions and 
using Eq. [O] The corrected results are shown for both orders 
of anisotropy in Fig.0] The error bars shown are statistical. The 
results for the first order anisotropy include contributions from 
a possible vertex shift. The upper and the lower limit of the 
boxes shown in Fig |4t show the asymmetric systematic error 
and correspond to no uncertainty and a maximum uncertainty 
in the position of the vertex, as discussed above. 

Fig-HJshows V2 as a function of centre of mass energy from 
various experiments (results mostly taken from compilation of 
Ref. 1 36 1). These have been measured in different experiments 
in different kinematic ranges using methods that vary in detail. 
The general behaviour shows a continuous increase in V2 as a 
function of centre of mass energy. The results of the present 



work are shown for two different centrality ranges correspond- 
ing to 10-20% of cross section and 10-30% of cross section. For 
comparison to the present results, values of V2 from PHOBOS, 
NA49 and CERES experiments |35 36 37 1 at nearly the same 
centre of mass energy are included. The various measurements 
are seen to agree well within errors. 



5 Anisotropy in Photons 

The photons incident on the PMD predominantly result from 
the two photon decay of the neutral pion, and if both pho- 
tons are detected in the PMD an additional apparent anisotropy 
will result from the kinematic correlation between the photons. 
The limited efficiency and purity of the detected photon sam- 
ple in the PMD affect the measured photon anisotropy values. 
The quantitative effect is estimated using simulations and is de- 
scribed in the following. 



5.1 Decay Effect 

The decay of neutral pions into two photons introduces correla- 
tions that can cause apparent anisotropics in the photon distri- 
butions which are greater than the actual anisotropy of the pi- 
ons. On the other hand, the process of decay smears the photon 
momenta relative to the initial pion momenta and can thereby 
dilute the initial correlation present in the neutral pions. The 
relative importance of these two competing effects has been 
shown to scale with the experimentally measured quantity Xm 
and enables a determination of the neutral pion anisotropy from 
the measured photon anisotropy 1 23 1 . However, the limited effi- 
ciency and the contamination of charged particles in the sample 
preclude a determination of the tt" anisotropy from the mea- 
sured 7 anisotropy in the present work. Using the anisotropy 
values of the charged particles folded with the tt" decay and 
experimental response allows one to determine the expected 
values of anisotropy of the photons in the PMD, as discussed 
below. 



5.2 Efficiency and Contamination 

The PMD records particle hits which include incident photons 
and a contamination of charged particle hits. These charged 
particles could be primary, or secondary rescattered particles. 
As noted above, the photon counting efficiency (e) of the sam- 
ple varies from 68% to 73% for central and peripheral events 
and the corresponding purity (p) of the sample varies from 65 % 
to 54%. 

The effect of decay, identification efficiency, and contam- 
ination on the observed photon anisotropy has been estimated 
using simulations. The simulations assume that the photon sam- 
ple contains a contribution from charged particle contamination 
which directly reflects the measured charged particle aniso- 
tropy in addition to a contribution from photons from 7r° de- 
cays with a tt" anisotropy which is also equal to the measured 
charged particle anisotropy. These simulations use the aniso- 
tropy values of the charged particles, vr*' decay kinematics, the 
PMD acceptance, and the purity of the PMD photon sample 
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to generate simulated data. The simulated sample is analyzed 
to obtain an estimate of the expected photon anisotropy corre- 
sponding to the observed charged particle anisotropy. 

The neutral pions were generated using the experimental 
pseudorapidity distribution of the charged pions 1 38 1 with an 
exponential p^ distribution (dN/dpy = p^ • exp(-6pT)). The 
7r° multiplicity values were chosen as half of those measured 
for charged particles with the SPMD for the same centrality 
selection. The second order anisotropy values of 7r° were cho- 
sen to be linearly increasing with px before saturating at p^ ^ 
1.5 GeV/c. For each centrality, the p^ dependence was chosen 
to reproduce the px integrated mean values of w„ of charged 
particles shown in Fig.|3 The linear dependence and the value 
of pt at saturation were both varied to estimate the system- 
atic errors. Neutral pions were generated and decayed and the 
decay photons were accepted if within the PMD acceptance. 
Using the measured photon multiplicity for a given centrality 
class, p ■ Nphotons photons were randomly selected from those 
falling onto the PMD, where p was assigned a value of 0.65 to 
0.54 corresponding to the centrality selection being simulated. 
A background contribution of ( 1 — p) • Nphotons charged parti- 
cles was added to the simulated event. This simulated data was 
then analyzed using the methods detailed in section 3.4. 

The systematic errors in the simulated results have been es- 
timated for both orders of anisotropy for each centrality. The 
contribution to the systematic error on the simulated results are 
shown for one particular centrality selection (20%-30% of min- 
imum bias cross section) in Table|2l 



Table 2. Contribution to systematic error from various sources for the 
simulated values of both orders of anisotropy of the 1^ photons distribu- 
tions for the centrality selection corresponding to 20-30% of minimum 
bias cross section. 



Source 


First Order 


Second order 




+0.005 


+0.002 


Charged Particle Anisotropy 








-0.004 


-0.003 


Purity of photon Sample 


±0.002 


±0.001 


Anisotropy of Contaminants 


-0.004 


-0.007 


rj and pr distribution of tt" 


± 0.001 


± 0.001 


Neutral Pion Multiplicity 


± 0.002 


± 0.001 



The relation between the neutral pion anisotropy and pho- 
ton anisotropy is not linear |23|. v^'^^^ has been estimated 
for a range of pion anisotropy values corresponding to the er- 
ror measured in this experiment. An increase (decrease) in the 
charged particle anisotropy increases (decreases) the anisotropy 
in the simulated results of v^'^^^ for all centralities. The per- 
centage change in w^*^^ corresponding to a 10% change in 
the purity of photon sample is small and is about the same for 
all centralities. If 50% of the contaminants are assumed to be 
isotropic, then the resulting anisotropy decreases for all cen- 
tralities, and contributes maximumally to the total systematic 
error The uncertainty in the neutral pion multiplicity also has 
a small effect, which is almost independent of centrality. 



The systematic error on the measured values of u^*^^ have 
been obtained by 

- increasing and decreasing the region of acceptance for the 
analysis 

- varying the size of the interval Aij between the two subevents 
in the range 0.03 to 0.07. 

- randomly removing up to 20% of the photons in the PMD. 

- obtaining the v^^^^ values by the correlation between the 
subevents as described in section 3.4 above. 

Repeating the analysis by rejecting clusters closer than twice 
the size of the scintillator pads did not change the w^*^^ val- 
ues. 



5.3 Results 

The measured values of w,^*^^ are shown in Fig.|6l The errors 
on w^*^^ are obtained by adding the systematic and the statis- 
tical error in quadrature and are also shown. The open triangles 
show the expected values of v^^^^ from the simulations de- 
scribed above. The statistical and systematic errors are added 
in quadrature and are shown by the shaded regions. The pho- 
ton anisotropy coefficients extracted from the simulated PMD 
data are consistent within errors with the measured PMD result. 
This demonstrates that the photon anisotropy results measured 
with the PMD are consistent with the charged particle results 
presented in Sec. 4. Note, however, that the results shown in 
Fig.6 include vr*' decay and charged particle contamination ef- 
fects, and should not be compared directly with other results. 
In particular, the use of the PMD photons themselves to de- 
termine the event plane gives a strong tt" decay effect. This is 
shown by the open circles in Fig.6 which show the simulation 
values of w^*^^ for the assumption of isotropic 7r° emission, 
Vn{TT)=0. It is seen that the first order photon anisotropy re- 
sults are dominated by the vr'^ decay effect while the second 
order photon anisotropy results also have a significant decay 
contribution, which cannot simply be subtracted from the mea- 
sured v|^*^^. Since the relationship between v^ and w^o de- 
pends in a non-trivial way on the 7r° (pT dependent) multiplic- 
ity and anisotropy |23 39 1 it is non-trivial to extract the 7r° 
anisotropy from the photon anisotropy measured with respect 
to the event plane determined using photons. It should be noted 
that the WA98 photon w„ results measured with the LEDA pho- 
ton detector covering the photon rapidity region 2.3 < y < 2.9 
II39I were obtained with respect to an event plane determined in 
the target fragmentation region. The v^*^^ measured using the 
PMD photons and w„ measured using the LEDA photons will 
therefore have different sensitivity to the tt" decay effect. The 
PMD results include effect of decay correlations, and therefore 
the current measurement, v^^^^, provides upper limits on the 
anisotropic flow coefficients u„. 



6 Summary 

The azimuthal angle distributions with respect to the event plane 
have been measured for charged particles and photons with 
full azimuthal coverage in the pseudorapidity region of 3.25 
< f] < 3.75 for 2"8pb + 208pb colHsions at 158 -A GeV/c. 
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A total of 0.25 million events, classified in seven centrality 
selections, have been used in the analysis. The Fourier coef- 
ficients of the azimuthal distributions have been extracted in 
several ways, all giving consistent results for the first order vi 
(directed) and V2 second order (elliptic) anisotropics for pho- 
tons and charged particles. The results show the expected trend 
of decreasing anisotropy with increasing centrality for both vi 
and V2 for charged particles and photons. Our results agree with 
the results reported by other experiments for near similar condi- 
tions. |35'36'37'39|. The observed anisotropics of the photon 
distributions compare well with those obtained from simula- 
tions that include the charged particle contamination and the 
correlations arising due to the decay of the neutral pions as- 
sumed to have the same anisotropics as measured for charged 
particles. 
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Fig. 1. Acceptance corrected distributions of (a) the first order event plane angle, i/>i, of the two subevents for the charged particles hits in the 
SPMD and the photons in the PMD in the region 3.25 < r] < 3.75 for the centrality class defined by 225.5 < Et < 298.6. (b)The same for 
second order event plane angle, ■ip'2. 
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Fig. 2. Distributions of azimuthal angles with respect to the first and second order event plane for the centrality class defined by 225.5 < Et 
< 298.6. (a) for charged particle hits in the SPMD and (b) for photons in the PMD, both in the pseudorapidity region 3.25 < rj < 3.75. 
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Fig. 3. The ratio v^^^" /vn for different number of occupied cells, x is the number of occupied cells scaled with the total number of cells. 
Different symbols correspond to different initial values of w„. The solid line represents Eg. II 31 
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Fig. 4. Anisotropy coefficients of the azimuthal distributions of charged particles in the pseudorapidity region 3.25 < r] < li.lS for different 
centralities characterised by the measured transverse energy, (a) First order, «i, where the shaded region indicates the extent of the total 
systematic error due to uncertainty in the vertex position, (b) second order, 112 . 



M.M. Aggarwal et al.: Azimuthal Anisotropy of Photon and Charged Particle Emission in ^"^^Pb + ^"^^Pb Collisions at 158 -A GeV/c 



13 



> 



0.06 
0.055 

0.05 
0.045 

0.04 
0.035 

0.03 
0.025 

0.02 
0.015 

0.01 



D 



WA98 charged particles(10%-30%) 
WA98 charged particles(10%-20%) 
CERES charged particles(21%-26%) 
NA49 pions(12.5%-33.5%) 



i 



+ 



^ 



10 



■I 



▲ AGS(E877) 

• STAR 

# Phenix 
ik Phobos 



10 



V^M (GeV) 



'NN 



Fig. 5. Second order coefficient for different centre of mass energies. The results of the present work are for charged particles and are shown 
for two different central! ty classes. The ordinates of CERES, NA49, STAR and PHENIX results have been shifted slightly for clarity. The error 
bars on WA98 points indicate the statistical and systematical errors added in quadrature. E877, STAR, PHENIX and PHOBOS V2 values are 
for charged particles. The centrality cuts for all experiments are comparable. 
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Fig. 6. (a) First order, Vi , and (b) second order, «2 , photon anisotropy coefficients in the pseudorapidity region 3.25 < r) < 3.75 for 
different centralities are shown by filled triangles. Statistical and systematical errors are added in quadrature and shown as bars on the filled 
triangles. Open triangles are the most probable values of v^^^ as expected from the simulation. The shaded regions indicate the simulation 
uncertainties as described in section 5.4. The open circles show the calculated values of u^'*^^ («„ (tt) — 0) assuming an isotropic distribution 
of pions with the dashed curve indicating a smooth polynomial fit to the open points. Note, however, that v^^^{v„{'iv) = 0) can not be 
directly subtracted from the v^'^'^{v„{it) > 0) to obtain the anisotropy flow coefficents u„, as explained in the text. 



